Large crystal local-field effects in the dynamical structure factor of rutile Ti02 
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We present ab initio time-dependent-density-functional calculations and non-resonant inelastic 
x-ray scattering measurements of the dynamical structure factor of rutile Ti02- Our calculations 
are in good agreement with experiment and prove the presence of large crystal local-field effects 
below the Ti M-edge, which yield a sharp loss peak at 14 eV whose intensity features a remarkable 
non-monotonic dependence on the wave vector. These effects, which impact the excitation spectra 
in the oxide more dramatically than in transition metals, provide a signature of the underlying 
electronic structure. 
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Titanium dioxide (Ti02) has been studied extensively 
for its remarkable electric, magnetic, catalytic, and elec- 
trochemical properties. Based on these properties, Ti02 
has been used in a wide variety of technological applica- 
tions, such as dielectric material for integrated electron- 
ics and photocatalyst for the decomposition of organic 
compounds^ 

In the last decade, many investigations have been car- 
ried out focusing on the structural, electronic, and op- 
tical properties of this transition-metal oxide. From the 
experimental point of view it has been studied by var- 
ious techniques, such as ultraviolet photoemissionf 2 * 3 . x- 
ray emission^ x-ray photoemission; 5 ^ Auger-electron^ 
and electron-energy-loss spectroscopies A^iiSi Theoretical 
investigations of this wide-band-gap semiconductor in- 
clude ab initio calculations of its structural, electronic, 
and optical propertiespii* 1 ^ 

Recently, Launay et a/. 10 have shown evidence of a 
rutile-phase characteristic peak in the low-energy loss 
spectra of Ti02, and Vast et a/— have investigated the 
impact of crystal local-field effects (LFE) on the electron 
energy-loss spectra of rutile Ti02- These authors used 
electron energy loss spectroscopy (EELS) to study the 
low wave vector behaviour of the energy-loss function of 
this material. They combined their experimental spectra 
with ab initio calculations, and concluded that local-field 
effects are only relevant at energies above 40 eV, where 
excitations from the Ti 3p semicore levels occur. 

In this paper, we report ab initio time-dependent- 
density-functional (TDDFT) calculations of the dynam- 
ical structure factor of rutile Ti02 at large momentum 
transfers. We also report non-resonant inelastic x-ray 
scattering (IXS) measurements, which are suited to study 
the short-wavelength behaviour of this quantity, and we 
observe large crystal local-field effects at low energies. 
These effects are absent at the small momentum trans- 
fers accesible by EELS; 13 however, we find that at larger 



momentum transfers, where inhomogeneities due to lo- 
calized d-states can be sampled, crystal local-field effects 
yield a sharp non-monotonic loss peak at 14 eV. 

Within the first Born approximation, the inelastic scat- 
tering cross-section for x-rays to transfer momentum 
7i(q + G) and energy Tilu to a periodic solid is charac- 
terized by the dynamical structure factor of the solid: 
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where (ei,e/) and (wi,uif) refer to the polarization vec- 
tor and frequency of the incident and scattered photon, 
respectively, the wave vector q is in the first Brillouin 
zone (BZ), and G is a vector of the reciprocal lattice. 
5(q + G,ui) represents the dynamical structure factor— 



5(q + G, lo) = -2 Q ImxG,G(q, w). 



(2) 



Here, XG,G'(q,w) is the density-response function of the 
solid, and f2 is the normalization volume. 

In the framework of TDDFT, 15 the exact density- 
response matrix of a periodic solid of electron density 
no(r) can be written as 

XG,G'(q,u>) =Xc,G'(q> w ) + Xg,gM> u ) 

Gi,G2 

x {»G 1 (q)+/G?,G a [no](q,w)}xG B ,G'(q,w). (3) 

Here, /q g , [n] (q, to) are the Fourier coefficients of the 
functional derivative of the time-dependent exchange- 
correlation (XC) potential of TDDFT, and Xg G'(q> w ) 
represents the density-response matrix of noninteract- 
ing Kohn-Sham (KS) electrons^ which is obtained from 
the knowledge of the eigcnfunctions and eigenvalues 
of the single-particle Kohn-Sham equation of density- 
functional theory (DFT)pi£*i£ We solve this equation self- 
consistently in the local-density approximation (LDA), 
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with use of the Perdew-Zunger parametrizatior. 19 of the 
Ceperley- Alder XC energy of a uniform electron gasS In 
the random-phase approximation (RPA), the XC kernel 
f^ c G , [n] (q, u>) entering Eq. © is taken to be zero; in the 
so-called adiabatic LDA (ALDA), it is approximated by 
an adiabatic local kernel of the form 
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where Vxc (n) is the XC potential of a homogeneous elec- 
tron gas of density n. 

In the calculations presented below, we first expand 
the LDA single-particle Bloch states in a linearized aug- 
mented plane wave (LAPW) basis,- 21 by dividing the unit 
cell into non-overlapping atomic spheres and the intersti- 
tial region between them. Inside the atomic spheres, the 
LAPW wave functions are expanded in spherical har- 
monics with Z max = 10. The Ti 3s and 3p states and 
the O 2p states are treated as semi-core states, and the 
Kohn-Sham equation is solved with a cutoff parameter 
RMTK max = 8 and an energy cutoff of 7.5 Ry. The 
Kohn-Sham density-response matrix has been computed 
as in Ref. I22I with a damping parameter of 0.2 eV. 

Non-resonant inelastic x-ray scattering measurements 
of the dynamical structure factor of rutile TiC>2 were ob- 
tained using the UNICAT undulator Beam Line on the 
Advanced Photon Source at Argonne National Labora- 
tory. The measurements were made in reflection geom- 
etry, with 7.6 keV x-rays, a range of wave vectors from 
0.5 to 3.98 A in the [001] direction, and an energy 
resolution of 1.1 eV. The surprising result of these mea- 
surements is the non-monotonic wave- vector dependence 
of the intensity of a sharp loss peak at 14 eV, which is 
absent in existing optical (q=0) and EELS measurements 
on this material and whose physics is qualitatively differ- 
ent from the well-known collective excitations in simple 
metals. Our ab initio calculations indicate that this fea- 
ture is the result of large crystal local-field effects near 
the plasma frequency. 

In order to investigate the impact of crystal local- 
field effects on the dynamical structure factor of rutile 
TiC>2, we have first neglected these effects by considering 
only the diagonal elements of the Kohn-Sham density- 
response matrix entering Eq. @ (diagonal calculation), 
and we have then solved this matrix equation with a given 
number of G vectors (full calculation) m& Well converged 
results have been obtained with the matrix size ranging 
from 63 x 63 to 113 x 113, depending on the momentum 
transfer and the energy range considered, and using 30 
points in the irreducible BZ (IBZ). 

The energy-loss function Im [— Cq (q, cj)p^ of rutile 
TiC>2 for vanishing momentum transfer (long wave- 
lengths) has been reported by Mo and Chingi^ with no 
inclusion of crystal local- field effects and more recently by 
Vast et alm& with full inclusion of these effects. At small 
wave vectors, two principal structures can be indentified 
in the energy-loss function at energies below the titanium 



M-edge: a plasmon peak at 12 eV, and a broad collective 
excitation at about 25 eV, stemming from the building up 
of collective modes of the strongly hybridized Ti 3c? and 
O 2p bands. Vast et alm& also reported loss-function cal- 
culations of rutile Ti02 for |q| w 0.4 A , showing that 
at these small wave vectors crystal local-field effects have 
little impact at energies below the titanium M-edge but 
drastically reduce the peak heights above the Ti M-edge. 

For vanishing and small momentum transfers our cal- 
culations esentially reproduce the calculations reported 
in Refs. [T3 and [I3 Besides, for large wave vectors not 
addressed before we find large crystal local-field effects, 
as shown by the dynamical structure factor per unit vol- 
ume s(q, lu) = S(q, u>)/Q of rutile TiC>2 which we plot in 
Fig. n for two values of the wave vector and for energies 
below the Ti M-edge. IXS measurements of this quan- 
tity are scaled to absolute units by first performing an 
absolute measurement on Al (with the use of the / sum- 
rule)^ and then scaling by the ratio of absorption coef- 
ficients for Al and Ti02- Considering the fact that the 
comparison between the calculations and measurements 
is in absolute units, with no adjustable parameters, the 
agreement between theory and experiment is remarkable. 

The top panels of Fig. show that in the absence 
of crystal local-field effects (dotted lines) a simple plas- 
mon peak is present at 12 eV, which is well defined for 

wave vectors up to a critical value of ~ 1.5 A where 
the plasmon excitation enters the continuum of intra- 
band particle- hole excitations. In the presence of local- 
field effects (solid lines), this plasmon is still present at 
|q| = 1.06 A (see Fig.^a)) but it is completely damped 
at larger wave vectors (see Fig.^b)); in addition, a sharp 
loss peak emerges at 14 eV, which brings our calculated 
energy-loss function into good agreement with experi- 
ment. 

We have also carried out calculations and measure- 
ments of the dynamical structure factor for larger val- 
ues of the momentum transfer. We have found that the 
new peak at 14 eV, which only appears at non-zero wave 
vectors, disappears at surprisingly large wave vectors of 
the order of 3 A without significant broadening. The 
oscillator strength of this peak, which is qualitatively dif- 
ferent from the collective excitatios in simple materials, 
thus features a remarkable, non-monotonic dependence 
on the wave vector. 

With the aim of establishing the nature of the large 
crystal local-field effects that yield the characteristic 
energy-loss function at 14 eV, we have considered the 
macroscopic dielectric function eA/(q, w), which within 
the first BZ is £M(q,w) = l/eg o(q, oj). Collective ex- 
citations occur at energies where both the real and the 
imaginary part of this dielectric function are close to zero. 

The bottom and central panels of Fig^exhibit the real 
and imaginary parts of eM(q 7 ^)- At the smallest wave 
vectors (not plotted here), crystal local- field effects are 
small and the real part of the macroscopic dielectric func- 
tion is zero at 12 eV where the imaginary part is small, 



3 




FIG. 1: Dynamical structure factor per unit volume s(q, u>) (top panel) and macroscopic dielectric function €m{<1,u) (middle 
and bottom panels) of Ti0 2 at (a) |q| = 1.06 A" 1 and (b) |q| = 1.59 A" 1 along the [001] direction. The solid (dotted) line 
represents the calculated ALDA spectrum with (without) LFE. The dashed line denotes the calculated RPA result with LFE. 
Open circles shows the IXS measurements normalized to the same absolute units as the theoretical spectra. 



thereby yielding a plasmon peak in the corresponding 
energy-loss function. However, as the wave vector in- 
creases, the rapid variation of microscopic electric fields 
acting on localized d states lead to a dramatic redistribu- 
tion of the strength in the imaginary part of the macro- 
scopic dielectric function (see middle panels of Fig. 0. As 
a result, collective excitations are considerably damped 
at 12 eV where Landau damping is now efficient. Instead, 
a well defined collective oscillation occurs just above the 
new continuum of particle- hole excitations. Indeed, the 
Kramers-Kronig relation between the real and imaginary 
parts of the macroscopic dielectric function guarantees 
that the real part has a dip whenever the imaginary part 
has a jump, and the result is a pronounced energy-loss 
peak (see upper panels of Fig. that is responsible for 
the realization of a surprisingly sharp collective mode 
at 14 eV. We have also carried out calculations of the 
dynamical structure factor of VO2 (in the rutile struc- 
ture), which also combines strong inhomogeneities of the 
Fermi sea with the presence of localized d states below 
and above the Fermi level, and we have found that this 
material also exhibits a non-monotonic sharp collective 



mode at low energies and large momentum transfers. 

We close this paper with a comparison of ALDA calcu- 
lations and IXS measurements of the dynamical structure 
factor of TiC>2 at energies both below and above the Ti 
M-edge. Figure [21 represents the results we have obtained 
for s(q, uj) with |q| = 1.32 A , showing that when crys- 
tal local-field effects are included the overall agreement 
between theory and experiment is very good. 

At low energies below the Ti M-edge, which is located 
at ~ 40 eV, crystal local-field effects yield a well-defined 
collective excitation at 14 eV, as discussed before. At 
energies above the Ti M-edge, a prominent structure is 
visible, which can be attributed to transitions from the 
occupied Ti semicore 3p states to the lowest conduction 
bands. In the absence of local-field effects (dotted line), 
the onset of these transitions would occur at 36 eV. How- 
ever, local-field effects bring the onset of the M-edge to 
~ 42 eV and reduces the corresponding peak heights, in 
close agreement with experiment. Nevertheless, there is 
still a small mismatch between the converged ALDA cal- 
culation and the measured spectrum at the position of 
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FIG. 2: Dynamical structure factor of Ti02 at |q| = 1.32 A 
along the [001] direction. The solid (dotted) line represents 
the calculated ALDA spectrum with (without) crystal local- 
field effects. The dashed line represents the calculated RPA 
spectrum with crystal local-field effects included. Open circles 
show the IXS measurements normalized to the same absolute 
units as the theoretical spectra. 



the main semicore peak. This mismatch stems from the 
fact that our calculated LDA ground state exhibits an 
upper semicore edge that is located a few eV higher than 
observed experimentally. Fig. [21 also shows that many- 
body XC effects not included in the RPA (dashed line) 
do not affect the low-energy collective excitation at 14 eV. 



These effects, however, slightly modify the structure of 
the collective excitation at ^25 eV (see also Figs, ^a) 
andQJb)) and yield a small shift of the semicore excita- 
tion onset and main peak towards lower frequencies. 

In summary, we have presented a combined theoretical 
and experimental investigation of the dynamical struc- 
ture factor of rutile TiC>2 at large momentum transfers. 
We have found that large crystal local-held effects in this 
wide-band-gap semiconductor yield a sharp loss peak at 
14 eV. An analysis of the macroscopic dielectric func- 
tion of this material leads us to the conclusion that this 
loss peak is originated in a collective excitation which re- 
mains present at surprisingly large momentum transfers 
and whose spectral weight features a remarkable non- 
monotonic dependence on the wave vector. This fea- 
ture is exhibited by both TDDFT calculations and IXS 
measurements, which show a remarkable agreement when 
presented in absolute units with no adjustable parame- 
ters. 
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